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Abstract 

Lattice radial quantization is introduced as a nonperturbative method intended to numeri- 
cally solve Euclidean conformal field theories that can be realized as fixed points of known 
Lagrangians. As an example, we employ a lattice shaped as a cylinder with a 2D Icosahedral 
cross-section to discretize dilatations in the 3D Ising model. Using this method, we obtain 
the preliminary estimate 77 = 0.034(10). 
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1 Introduction 



Conformal field theories are theoretically interesting in their own right, have applications in 
Condensed Matter and Statistical Physics and are relevant to Particle Physics in general. They 
may play a central role in the yet unknown correct description of Particle Physics beyond the 
Standard Model. Traditional numerical methods work well when the correlation length can be 
kept at a few lattice spacings but are ill suited for conformal field theories. Here we propose an 
untraditional numerical formulation, based on radial quantization. 

A classical local Euclidean action on with no dimensional parameters can be rewrit- 
ten in polar coordinates with radius r as a field theory on S'^^^ x R with invariance under shifts 
along the factor M parametrized by log(r) [l]. The shift symmetry corresponds to dilatations 
about the chosen origin in R"^. The corresponding generator is the dilatation operator and its 
eigenvalues are operator dimensions. The critical new feature of our proposal is to introduce 
a uniform lattice discretization of log(r); this provides a reasonable match between numerical 
resources and the relevant degrees of freedom. A lattice of linear size L in log(r) in our proposal 
replaces a lattice of linear size when R*^ is regularized directly. Using the 3D Ising model, we 
carry out a modest numerical test of lattice radial quantization of the Wilson-Fisher fixed point 
field theory. 



2 Preserving dilatation invariance 

It may sound surprising that one can regularize a field theory that is classically conformally 
invariant while preserving dilatations because this seems to say that an asymptotically free 
theory, such as QCD, does not necessarily generate a mass. We first clear up this point by 
studying the large N limit of the 0{N) sigma model in 2D, and in the course of doing so we 
learn how to identify a theory that does stay conformal after quantization. 

2.1 An example with asymptotic freedom 

Consider the quantization of the two dimensional 0{N)/0{N — 1) nonlinear sigma model at 
large iV, where it can be solved exactly. One can choose to preserve either dilatations together 
with rotations or translations in two directions, but not all of these symmetries. If one forces 
the preservation of all these symmetries in the quantum theory, it would end up being free. The 
partition function of the model is 

Z = j [d0][da]e-^ , S = ^jcfx [{da(t>? + i<T[<t}^ - N/X]] , (1) 

where cf) = {(pi, (/)2, - ■ ■ jfpN) is the N-component vector field and the limit iV — t- oo is taken at 
fixed A, with A = g'^N. The auxiliary field a enforces the constraint cfp' = (j) ■ (f) = l/g^ . 
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2.2 Standard, translation invariant analysis 



At large A^, Z can be computed by saddle point. An implicit cutoff is assumed - it will be made 
explicit later on. In the standard approach, the global 0{N) symmetry, which is nonlinearly 
realized at the classical level, becomes linearly realized at the quantum level by a multiplet of 
N massive scalar particles. Their mass, is obtained by seeking a translation invariant saddle 
point of the form cr{x) = —ifjp. In Fourier space we have at the saddle: 

0(^) = / d'pe''"'4>{p) ^S=\j ^0 + ^^)~^^p) _ (2) 

with ip) = (f){—p). After integrating over the scalar fields we obtain Z = J diJ?e~^^^^^^'^\ The 
value of /i^ is obtained by minimizing F{ii'^, ^) — \ I ln(p^ + fj?) — jxfJ-'^ and is /i = ^ ^^^^ ^ 

for a UV cutoff of the form < A^, where A is a large momentum cutoff of mass dimension. 
For large A/fi we have A ~ in(A/fj.) • This is the standard result with dimensional A and /i. 



2.3 Dilatation invariant analysis 



Following [T] we first go to radial coordinates, xi = rcosO, X2 = rshi9 and then to Fourier 
space in 6: 



■m=— oo 



and 



^ oo 
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r e 
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-rn{r) 



27r ^ 

m=— oo 



, 0-m{r) = (T*^{r) . 



(3) 



(4) 



We now map the radial coordinate to the infinite line. This step turns dilatations into transla- 
tions in t and inversion into a parity in t. 



r = e*, dr = e*dt, ^m{t) = cl)„^{r), ^^(t) = f7m(r) 



(5) 



In the new variables the classical action is 



1 

in 



2A 



-N 



dtSo(t) 



(6) 



with 



-^+"1 j dm+m',0 + -^^~m,-m'{t) 



6it - t') 



(7) 



We now look for a rotation invariant saddle point that also is translation invariant in t, = 
27rl](5m,o- In Cartesian coordinates, such a saddle is scale invariant, but singles out one point, 
serving as the origin. 
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After integrating out the vector field, we obtain the following unregulated saddle point 



equation: 
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The sum over m is unrestricted and can be done in closed form, producing 
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(8) 



(9) 



We now set iS = /i^ and choose a cutoff A > 0. A > is dimensionless since t was dimensionless 
and k is the conjugate momentum to t. We regulate by integrating over — A < < A assuming 
A//i, A//io ^ 1- This preserves also the t ^ —t symmetry. 
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(10) 



A/}, is a renormalized coupling, /i and are dimensionless. /^o is an arbitrary renormalization 
point. The bare coupling depends on the cutoff in the standard way, except the cutoff is dimen- 



sionless: A 



2tt 



A nonzero saddle point will always be found. The saddle point value ^ 



\n{2h/^JLo)■ 

never exceeds a maximal value, which depends on and is defined by 



^max(/^o) '■ 



dk 



coth vr \/^^^+l^^^^Jpo) 



(11) 



The propagator of the scalar field is 

dk 



Jk{t-t') 



27r k"^ + + /U^ 



I e 



(12) 



2.4 Cartesian Translation Invariance 



Reading off the set of dilatation eigenvalues displayed by eq. (fT^ we see that it does not contain 
an equally spaced ladder. Consequently, we cannot construct translation generators satisfying 
the correct commutation relations [2] with dilatations in the sector generated by the action of $ 
on the vacuum. The deviation of the dilatation spectrum from equal spacing is small if m ^ /i. 
There is approximate translation invariance for angular separations that are small relatively 
to -. Such separations can be resolved only by waves with high enough angular momenta m. 
Because inversion has also been preserved in the quantization, if translations could be realized, 
special conformal transformations would come in automatically and the full conformal group 
would be realized. Because of rotation invariance, only one linear combination of translations 
needs to be considered in detail. 
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Comparing the manipulations we did here to those we would have done had we de- 
clared we are interested in the traditional view of the model at finite non-zero temperature, 
we realize that while the language has changed, the equations have not. There exists only one 
bosonic vector 0{N)/0{N — 1) sigma model on § x M. In traditional terms, the circumference of 
S in units of the dynamically generated mass at zero temperature is a free parameter. Insisting 
that S is a geometric circle made its length equal to 27r in the untraditional view. In tradi- 
tional terms the propagator decays exponentially in the Cartesian spacetime directions, while 
in untraditional terms the propagator decays exponentially in the coordinate conjugate to the 
dilatation generator. 

2.5 One useful lesson 

We did not get anything truly new in the untraditional quantization; the "novelty" boiled down 
to semantics. But, we learned a lesson which will be useful in the following: if we started from a 
classical conformal theory that allowed this symmetry to be fully preserved at the quantum level 
and chose to preserve dilatations and rotations explicitly by the cutoff procedure, translations 
would be restored in the continuum limit. The generators of translations could be reconstructed 
from the commutation relations of the conformal algebra using the available generators of ro- 
tations, dilatation and inversion [2]. Eventually we would get the same theory we would have 
gotten from a traditional lattice regularization that preserved a discrete version of translations 
from the start but left dilatations and rotations unpreserved. However, at the lattice regulated 
level, the choice to preserve dilatations seems more efficient. 



3 Lattice radial quantization in three dimensions 

In the continuum we would be dealing with a theory on §^ x M. The critical exponents of the same 
CFT on can be extracted from the dilatation operator on S'^ x M [3] . In the two dimensional 
case one has § x M and the same can be done [1] . This works well because the factor § is easy to 
discretize. In three dimensions we need to discretize §^ and it becomes impossible to replace the 
continuum 0(3) symmetry by symmetries under ever larger discrete subgroups as the lattice is 
refined. One is limited to a few discrete subgroups. Hopefully, this gets compensated by their 
stronger, nonabelian structure. One can latticize S'^ [3] by the set of vertices of a platonic solid, 
but this does not work well [5]. A way to improve on this is to use regular lattice refinements 
of the faces of the platonic solid [6]. For a cube, this approach produced reasonably accurate 
exponents in the Z2 even sector of the Ising model. However, the estimate of the exponent 77 
in the odd sector had a substantially larger deviation from the known value [7]. Our numerical 
experiment will exploit the lesson we learned in the previous section on a lattice similar to that 
of [6]. 

To have as much symmetry as possible we employ a closed two dimensional surface 
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S, consisting of the boundary of an icosahedral solid instead of a cube. The hnear distances 
between points on S and the origin are bounded to a relatively narrow interval on the positive 
axis. The icosahedron has 12 vertices and 20 faces given by identical flat equilateral triangles. 
Its symmetry group 1^ is a 120 element subgroup of 0(3). The angular momenta I = 0,1,2 
representations of 0(3) remain irreducible representations under Ih. This is enough to isolate a 
scalar primary state and two of its immediate descendant states on the basis of their behavior 
under 1^. is visualized as a set of 20 infinite triangular pyramids glued together with the 
apexes meeting at the origin. In isolation, each pyramid can be (non-conformally) mapped 
into a uniform prism with an equilateral triangle as cross section. These prisms are then glued 
together in a manner inherited from the pyramids. We end up with S x M. 

We regularize S by replacing each face by a triangular piece of a regular two dimen- 
sional triangular lattice. Consequently, each side of the original triangle has extra s — 1 equally 
spaced points placed on it (s > 1). The number of sites on one icosahedral shell is 2 + lOs^. 
Identical grids replacing S are strung along an equally spaced infinite straight line, making up 
a cylinder. Matching vertices on neighboring S-slices are connected by links of equal lengths. 

For a finite cylinder, the partition function is of the usual form, 

where {ij) denotes a nearest-neighbor pair on the same icosahedral shell and t = 0, ■ ■ ■ T — 1 
labels the shells along the cylinder. The trace is a sum over the Ising spin values cr{t,i) = ±1 
for all sites (t, i). 

To get to the continuum limit of the Wilson-Fisher conformal field theory, we need to 
tune /3 to its infinite s critical point. The coupling /3 has been chosen the same in all directions 
for simplicity. The intra-shell lattice spacing and the lattice spacing on the triangulated shells 
are fixed by demanding that descendants have equal unit spacing as we approach the continuum 
limit. This is a self-consistent definition only if the descendant "masses" fall on a straight line. 
This happens only at the critical value of /3. To get information on the spectrum of the transfer 
matrix it is convenient to compactify the infinite axis of the cylinder to a circle with periodic 
boundary conditions on the spins. We chose our cylinders to have lengths T which scale with 
the refinement as T = ps. 

Our focus is on the 3D Ising critical exponent rj. This number is determined by the 
dimension A^- of the lowest primary in the Z2 odd sector of the space the transfer matrix acts on. 
To get the primary state and some of its descendants we define field variables on each t-labeled 
shell: 

a{t, l,m) = J2 ^^rYriOu 4>iMt, ei,(t>,) . (14) 

i 

S is parameterized by the polar and azimuthal angles 9, (j). t corresponds to the logarithm of 
the amount of dilatation a specific S embedded in would have undergone. The Y^{9, (j)) are 
the standard spherical harmonics with / = 0, 1, 2. The measure Afi, is computed from the area 
of the spherical triangles adjacent to 9i,(j)i projecting onto a unit sphere in order to assure a 
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rapid convergence to the orthonormality condition for the yj™'(0, (/>)'s in the continuum. The Ih 
symmetries and dilatations guarantee (ait, I, m)a{t', /', m')) = Ci{t — t')6ii'6mm' for /, /' = 0, 1, 2 
with \m\ < I , \m'\ <l'. From the / = channel we get the dimension of the primary state and 
from the I = 1,2-channels the dimensions of the next two descendant states. We denote these 
three eigenvalues of the dilatation operator by A/, where A/ plays the role of a mass when t is 
regarded as a time variable. The eigenvalues of the lattice dilatation operator corresponding to 
the continuum dimensionless numbers A; are denoted by /ij. The theory now has an ultraviolet 
cutoff denoted by A. is dimensionless, determining the smallest angular separation on the 
discretized S. 

A, = (15) 

A oc s as s — )• oo and we anticipate an integer spacing of descendants, A/ = Ao + ^, in the s — )• oo 
limit. Having three values of at our disposal due to the icosahedral symmetry, we can test 
numerically whether they indeed fall on one straight line. The slope of this line determines A 
and fixes the ratio of the lattice spacing on the sphere relative to radial axis in t. The will go 
as at large s and the standard critical exponent rj is given by the anomalous portion of A/ 
for any h r] = 2{Ai - I - 1/2). 

It is known that r] is small relative to 1 [7] • This puts a high accuracy demand on the 
Hi^s, which are the directly measurable quantities. 



4 Results 



First the critical point /3c was determined by computing the 4th-order Binder cumulant [8], 

UiP,^,P) = -^-j^- (16) 

Our initial computations of moments of the magnetization were performed using the Swendsen- 
Wang cluster algorithm [9]. Later, we switched to the more efficient [10] single cluster Wolff 
algorithm [TT]. To locate the critical point, we used aspect ratios p = 4,8. with s > 16 and /3 
such that 1/3 - ^ds'^/" < 0.0012. We then fit the data to: 

U{(3,s,p) ~ C/(/3c,oo,p) +ai(p)[/3 - Pds'/" + h{p)s-'^. (17) 

Subsets of the data used appear in Fig. [H The best fit gave /3c = 0.16098691(1) with critical 
Binder cumulants C/(/3c,oo,4) = 0.3032(1) and C/(/3c,oo,8) = 0.18765(5) with a xVdof = 1.37, 
dof = 212. The best-fit exponents v = 0.621(4) and u = 0.797(6) are consistent with the values 
in [7] at the 2-sigma level. 

To compute the correlation functions, C/(i) = {a{t + tQ,l,m)a{tQ,l,m)) , we generated 
ensembles at /3 = 0.160987. Our final results for /3c from the Binder cumulant suggest this 
might be slightly in the ordered phase as s — t- oo. Each independent run was thermalized using 
2048 sweeps of the Wolff algorithm followed by 8192 sweeps with one estimate of the correlation 
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Figure 1: Determining /5c from the Binder cumulants U{s) = 1 — (M^)/3(M^)^ near the pseudo- 
critical point for two different aspect ratios p and increasing values of s. 



function after each sweep. We defined a sweep to be 19s/2 Wolff cluster updates. This choice 
sets the average number of spins flipped per sweep to be about equal to the total volume. All 
results for a given run are then averaged together to form a single blocked estimate and many 
thousands of independent block estimates are combined to form the ensemble at each s, as shown 
in Table [TJ The jackknife method was used to estimate errors. 
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Table 1: is the number of independent runs used for the correlation function at each refinement 
s. Other details of the runs are described in the text. 

We found it very useful to evaluate the correlation function using the momentum- 
space single cluster improved estimator method [12\ I13j . As we are only interested in the 
rotation-invariant part of the correlation function on any given lattice, we also average over the 
m = states in the various representations labeled by /. Since our value of /3 should 

be slightly larger than the pseudo-critical coupling at any finite s, we expect that our I = 
correlation function will have a small disconnected contribution, which should vanish as s — )• cx), 
assuming /? = /3c. If the disconnected contribution remained non-zero, it would suggest that 
P > 13c] otherwise, we would conclude that 13 < /3c- The contribution of the lowest eigenstate of 
the transfer matrix to the connected correlation function in momentum space is given by 



Ci{k) = co5ifi5k,o + ai 



(l - e-'^'^) sinh(^i) 
sinh2(/ii/2) -hsin2(A;/2) 



(18) 
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where k = Inn/T with n = 0, ■ ■ ■ T — 1 is the momentum conjugate to t along the cylinder 
and we have included the disconnected contribution cq. We can determine cq by fitting Co (A;) 
for k ^ and subtracting a smooth extrapolation of liiRk^Q C(){k). We found that our data 
required parameterizing the lowest eigenvalue plus at least two higher eigenvalues. We got 
excellent fits with x^/dof < 1 and estimates of the lowest eigenvalues which are essentially free 
of contamination. The result is shown in Fig. [2j A linear extrapolation in 1 /s indicates that 
the disconnected contribution vanishes at finite s ~ 64. Then, /3 = 0.160987 < /3c, in apparent 
contradiction to the Binder cumulant result. 

Icosahedron(s), T = 8xs, (3 = 0.16098700 

2x10"^ 
1x10"^ 
o° 0x10"^ 
-1x10""^ 
-2x10"^ 

0.00 0.05 0.10 

1/s 

Figure 2: The disconnected contribution to the correlation function as described in the text. 

We relate the /U;'s for I = 0,1,2 to the eigenvalues of the dilatation operator by 
m = A~^[Ao + /] where ~ ci/s as s — t- oo. We find ci ~ 1.51(1) with the uncertainty 
dominated by systematic error. The the left figure of Fig. [3] shows evidence for also sub-leading 
contributions, 0{l/s'^). We confirm the equal spacing rule of descendants by examining the 
ratios, {^1+2 — fJ-i+i) / ifJ'i+i — as shown on the right in Fig. [3l We are now justified to 
estimate the scaling dimension of the primary operator using 

Ao = 

2 

as shown in Fig. [H The determination is independent of A. 

A quite conservative estimate is r] = 0.034(10), consistent with other estimates [7]. 
Our numerical study could be extended to include additional primary operators in both the Z2- 
odd and Z2-even sectors as well as a direct test of the restoration of full conformal symmetry for 
2- and 3-point correlators. However as a cautionary note, we now briefiy describe the behavior 
of the eigenvalues for the descendants with / = 3. 

The 7-dimensional / = 3 irrep of 0(3) splits under Ih into the sum of the 4-dimensional 
irrep G and the 3-dimensional irrep T2. The / = 3 states m = ±2, ±3 appear in both irreps of 




fii + Hi/ 1 + I' 



l-V 



(19) 
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Figure 3: The left figure shows the scahng of A with s. The extrapolated value is 1.51(1) with 
the uncertainty dominated by the systematic difference between the two estimates. The right 
figure checks whether the m fall on a straight line. We fit to a linear function and find an 
intercept of 0.994(4) and slope of 0.0(2) with xVdof = 0.43 for 11 dof. 



Ih but the states m = ±1 only appear in G, while the state m = only appears in T2. In Fig. [5] 
we show the / = 3 eigenvalues labeled by m = 0, ±1. If they extrapolated to the same value as 
s — >• 00, it would be an indication that 0(3) symmetry is restored in the continuum limit. The 
numerical evidence so far suggests the contrary, indicating that the continuum limit with the 
icosahedral shell is not fully reproducing the conformal fixed Wilson-Fisher fixed point of the 
3D Ising model on flat three-space. More reflned methods are needed, as suggested below. 



5 Discussion and Outlook 



The method of lattice radial quantization presented in this paper holds promise as a practical 
nonperturbative tool for conformal field theory. We plan additional numerical and theoretical 
studies to realize this goal. 

An important theoretical question remains: How do exponents associated with / = 
0, 1, 2 in our icosahedral shell system at s — )• 00 relate in principle to their exact values in the 
3D Ising model? The numbers in [6] and our ?7-value are consistent with equality, but are far 
from providing overwhelming numerical evidence. Indeed, the apparent lack of convergence to a 
single 0(3) irreducible representation for / = 3 in Fig. [5] suggests that the continuum limit of our 
radial-icosahedral model is not equivalent to the 3D Ising model at the Wilson-Fisher fixed point 
on flat three space. A failure to reproduce the correct spectral degeneracy o{ I > 2 descendants 
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Figure 4: The scaling exponent of lowest Z2-odd primary operator vs. 1/s. The lower (red) 
points come from //q and fii. The determination of fiQ is unreliable because of the disconnected 
contribution to the correlator, given the uncertainty whether /3c < 0.160987 Pc > 0.160987. A 
linear extrapolation from upper (blue) points, based on fj,i and /X2, gives Aq = 0.517(5), which 
is consistent with the estimate 0.5182(3) in [7]. 



at s — )■ cxD would indicate distorted spectra for the primaries as well. The icosahedral shell may 
be sufficiently spherical to keep these distortions in r/ to a level below our current estimate. 

We are now studying in the continuum the consequence of quantization on S x M where 
S is sectionally flat. Also, we are implementing improved actions for our lattice representation 
with nearest neighbor bonds weighted to approach the continuum metric on x M, similar to 
the discretization employed for our single spin operators. Both theoretical considerations and 
numerical tests are needed to clarify how a radial lattice can reproduce a conformal Euclidean 
field theory on flat space in the continuum limit. 

Lattice radial quantization should apply also to four dimensional gauge theories with 
an amount of matter where the IR behavior is suspected to be controlled by an interacting 
conformal theory. Historically, interest in such theories was motivated by a search for realizations 
of walking Technicolor scenarios |14j in their vicinity. Independently of the relevance of the latter 
to Nature, there is theoretical interest in these theories, in particular in the large N limit for 
the gauge group SU{N) because of the potential existence of a specific string dual in the IR. In 
the case of 4D gauge theories, the anomalous dimensions of interest are typically large, so the 
demands on numerical accuracy are weaker than those for r]. 



Acknowledgments: RGB acknowledges support under DOE grants DE-FG02-91ER40676, 
DE-FC02-06ER41440 and NSF grants OCI-0749317, OCI-0749202. RGB has benefited from 



11 



Icosahedron(s), T = 8xs, (3 = 0.16098700 




1.60 I ' ' 1 ' ' ' ' 1 ' ' 1 

0.00 0,05 0.10 

1/s 
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difference is a measure of the breaking of 0(3) symmetry at that s. 
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